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Abstract: 

Using the generalized Gutzwiller method we present results on the ferromagnetic 
behavior of extended Hubbard models with two degenerate d(e g ) orbitals. We 
find significant differences to results obtained from Hartree-Fock theory. 
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Recently, the Gutzwiller variational method has been generalized for n > 2 
correlated orbitals per site In this paper we present studies of the itinerant 
ferromagnetism using a two-band Hubbard model with degenerate orbitals and 
general on-site interactions. 

The one-band Gutzwiller variational wavefunction is given by ||] 

|*> = g d |* > = II I 1 " U " 9)D S ] \^ ) , (1) 

8=1 

where l^o) is an uncorrelated wavefunction on a lattice of L sites and the oper- 
ators D s = n s ^n s [ measure the double occupancy of sites s. By combinatorics, 
one is led to the hopping 'loss factors' 

1 r ]2 

^ = n \ v m k m o + Vm k imi , (2) 



m k (l - m k ) 

where m , mi, m 2 , m 12 (mo, mi, m"2,mi2) are the respective gross (net) occu- 
pancies of sites: empty, spin-up, spin-down, double. There are relations between 
gross and net occupancies such as mi = mi— mi2 or mo = m"o = 1— mi— m 2 — rfi\2- 

The extension to cases of n > 2 correlated orbitals (a, 6, ..) leads to J n = 4 n 
different possibilities for a single site occupancy, and K n = J n — (2n + 1) of them 
represent multiple ones. The basic idea of our generalization was to include all 
K n multiple occupancy operators. 

Using the notation af — 1, aj=2, b]= 3, 6|— 4, the resulting generalized 
hopping loss factors have been found to be 

qkk = — jzr r[Vm k fn + '^m kl mi (3) 



+ Yl v mkl P mi p + Y1 v m ki Pq m P q + • • •] 2 

9fci = 7m9u • (4) 

It is now possible to investigate extensions of the Hubbard model for arbitrary 



numbers of orbitals a, (3 and more general on-site interactions: 



H 




(5) 



For our present study we have focussed on a ferromagnetic case; i.e. we have 
chosen iV-particle wave functions , |\l/o) which allow that mi = 1713 = m + = 
(n-j-) 7^ (ni) = m_. This spin-splitting is controlled by an additional variational 
parameter A (m + = m + A; m_ = m — A), representing the magnetization. We 
have studied the n = 2 Hamiltonian with two e g — type orbitals on a simple cubic 
lattice and have included first nearest neighbor (INN) and 2NN hopping terms 



given by §: T dda (lNN) = leV, T dda (2NN) = 0.25eV, and T dd6 : T ddn : T dda = 



0.1 : —0.3 : 1. Further, there are three interaction parameters U aa = U, U a @ = 
U' and J a/3 = J. We employ the relation 2J = U — U', which is valid for e g — 
orbitals in the limit of vanishing configuration interaction [[|. 

There are seven variational parameters representing the multiple occupancies: 
/ = m 123 4, t+ = m 123 = m lu , t- = m 124 = m 2U , dt+ = mi3, d- = m 24 d d = 
mi2 = W34, g? = mi4 = m23. Further, there are two different loss factors 
Q+ = Qu = <?33 = 913 = Q31 and g_ = q 22 = qu = <?24 = <?42, with 





+ v / ^(v / / + V / ^))] 2 - 



The ground state energy function is 



E = q + e(m + ) + g_e(m_) 



(7) 



+2Ud d + 2U'd + ([/' - J)(d t+ + dt.) + 



+(2U + 4U' -2J)(t+ + *_ + /) 
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with e(m) being the kinetic energy of the uncorrelated bands. Note that the 
number of electrons per atom is given by e/a = 4m. We have studied band 
fillings, which, for the paramagnetic case, are close to the biggest peak in the 
single particle DOS. 

Fig. 1 shows the magnetization A as a function of U (with J = 0.2 U) for 
the values m = 0.35 (Fig la) and m = 0.3 (Fig. lb). Also shown are the values 
of N(E F ). When e/a is chosen so that N(E F ) of the uncorrelated case coincides 
with the DOS peak (m = 0.3), there exists a kind of Stoner criterion for the 
ferromagnetic instability. Then, A starts as a continuous function of (U, U', J) 
(Fig. lb). Away from this peak, there is a first order transition to a state of 
finite magnetization before the Stoner criterion is met (Fig. la). We observe 
further discontinuous increases of the magnetization with increasing interaction 
(U, U', J). These jumps are related to structures in the DOS which pass through 
the Fermi level with increasing splitting of majority and minority bands (see Fig. 
la). 

A Hartree-Fock treatment of the model leads to onsets of ferromagnetism at 
much smaller values of the interaction parameters U,U',J (Fig. 2). The first 
order transition from A = to finite A is also seen in HF theory; however no 
further jumps in magnetization are observed. 

The magnetic phase diagram (Fig. 3) elucidates the role of the on-site ex- 
change for the formation of itinerant ferromagnetism. 

In conclusion, our study of itinerant ferromagnetism indicates significant dif- 
ferences between the Gutzwiller method and the Hartree-Fock theory. We hope 
that the generalized Gutzwiller method will enable us to judge the quality of effec- 
tive one-particle theories for itinerant magnetism such as spin density functional 
theory for 3d transition metals. 
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Figures 



Figure 1: Fermi density of states N E f and interaction U (with J = 0.2 U) versus 
magnetization A for two bandfillings m. The dotted vertical lines in (a) indicate 
first order changes in A. 

Figure 2: Comparison of Gutzwiller and Hartree-Fock results for A versus U 
(J = 0.2 U). Dashed lines: m = 0.3, solid lines: m = 0.35. Inset shows the jump 
in A near U = 7.95 in the Gutzwiller result for m = 0.35. 

Figure 3: Phase diagram in the J/U versus U plane for the Gutzwiller results 
with m = 0.3 (note that U' = U-2J). 
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